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Abstract 

In this paper we review the theory of sihcon nanowires. We focus on nanowires 
with diameters below 10 nm, where quantum effects become important and the 
properties diverge significantly from those of bulk silicon. These wires can be effi- 
ciently treated within electronic structure simulation methods and will be among 
the most important functional blocks of future nanoelectronic devices. Firstly, we 
review the structural properties of silicon nanowires, emphasizing the close con- 
nection between the growth orientation, the cross-section and the bounding facets. 
Secondly, we discuss the electronic structure of pristine and doped nanowires, which 
hold the ultimate key for their applicability in novel electronic devices. Finally, we 
review transport properties where some of the most important limitations in the 
performances of nanowire-based devices can lay. Many of the unique properties of 
these systems are at the same time defying challenges and opportunities for great 
technological advances. 
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I. INTRODUCTION 



One-dimensional nanostructured systems have attracted a great attention in the last two 
decades, with this interest extraordinary boost ed by the faci le synthesis of carbon nanotubes 
(CNTs) reported in the beginning of the 90s ( lijimj . 1991 ). The reason is twofold: on the 
one hand they have proved to be an excellent test-bed to study the most intriguing physical 
effects, whereas on the other hand they are believed to be among the most important building 
blocks of the next generation of electronic devices. 

CNTs are hollow cylinders o btained by rolling up one or more graphene sheets, a one- 



atom-thick allotrope o: 



carbon ( 



Charlier et al. 



20071 ). The symmetry and the electronic 



structure of graphene (INeto et all l2009l ) are such that the properties of the CNT depends 
critically on the exact way it is wrapped up, and it can be either metallic or semiconducting. 
This confer the CNTs with a richer physics, but it is clearly far from ideal from the viewpoint 
of applications, especially when -as it is the case- a simple route to selectively grow one 
type of CNT or the other is lacking. 

Nanowires are an extremely attractive alternative to CNTs, because it is much eas- 
ier to control their electrical properties and, as long as the surface is properly passivated 
-something that occurs naturally during or right after growth-, they are invariably semi- 
conducting ^ 

Silicon nanowires (SiNWs), in particular, look like a very appealing choice, since they 
provide the ideal interface with the existing Si devices, while taking advantage from a 



tractable material techno^ 



nique (IWagner and Ellis 



ogy. SiNWs are comraonly 



1964 



grown by the vapor-liquid-solid tech- 



Westwater et all . 119971 ). where a Au nanoparticle is used to 



catalyze SiH4 decomposition. Briefly, the Au particle is deposited onto a Si surface and react 
with the Si atoms of the substrate, forming Au-Si alloy droplets. These droplets adsorb Si 



^ Below we will discuss explicitly a few cases where nanowires derived from semiconducting solids can be 
metallic, whereas nanowires made of metal atoms are beyond the scope of this paper and are not being 
discussed. 
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from the vapor phase, resulting in a supersaturated state where the Si atoms precipitate 
and the SiNW starts nucleating ^. 



As David K. Ferry illustrates in an enlightening paper (IFerryl . l2008l ). nanowires could 
provide the paradigm shift needed to continue improving the density and the performances 



of electronic circuits. For almost four decades t 



described by the well-known Moore's law (IMoore 



l e incr ease in computing power has been 



19651 ). which has been standing on three 



pillars: (a) the increase of the size of the microchips; (b) the reduction of the transistor size, 
and (c) the circuit cleverness, that is the reduction of the number of devices required to 
perform a certain function. While the first of these driving forces played a significant role 
only in the pioneering years of solid state electronics, the reduction of device size has a pivotal 
role, since the physical limit of material scaling is rapidly approaching. Nanowires can lead to 
an obvious benefit concerning the miniaturization, thanks to bottom-up growth that allows 
overcoming the limit of conventional lithography-based top-down design. Subtler are the 
perspective advantages concerning circuit cleverness, which can be significantly improved by 
taking advantage of the coexisting nature of interconnection and active device of nanowires. 
In particular, a replacement of metallic vias with vertical transistors is envisaged. These 
new circuits could be easily reconfigurat ed to perform different operations, achieving a much 



higher level of integration ([Ferry 



20081 ). Additionally, compared to classical planar device 



technology, nanowires can better accommodat e all-around ^aies (see Fig. [H), which improve 



field-effect efficiency and device performances (iColingd . 120041 : iNg et all 1200J) and mobilities 
of ~ 1000 cm^ V~^s~ ^ , substantially larger than those obtained in conventional Si devices, 
have been obtained ( ICui et all l2003l ) ^. 



Several pr omising app 



tron devices (jChung et al 



2006 



Hu et al. 



2008 



20041 ). logi c gates (jHuang et al. 



i cations have already been demonstrated, 



2000 



Lu et al. 



Cui and Liebei 



2008b 



tovoltaics (jKempa et al. 



2008 



2001 



Wang et al. 



Cui et al. 



2006b 



r anging from elec 



2003 



Yu et al 



200lh. non-vo 



Tian et al 



2009, 



atile m emories (IDi 



Goldberger et al. 



200ol: 



20071 ). photonics (iGudiksen et al. 



an et al. 



Zheng et al 



2002h . Dho; 



2002 : 



See 



Wang et al\ (|2008f ) for a comprehensive review of the growth techniques. 



It is difBcult to make rigorous comparisons, because the mobihty has a strong inverse dependence on 
the dopant density which is seldom known with accuracy i n nano wires. However, the peak value of 



Cui et al. 



(|2003l ) must be compared with the typical 



1350 cm^V^^s^^ obtained for the hole mobility by 
values for bulk Si of ^ 400 cni^V^^s^^ and ^ 100 cm^V^^s^^ for an acceptor concentration of 10^^ cm~'^ 
and 10^^ cm~^, respectively. 
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Pauzauskie and Yand. 120061) . to biological sensors (ICui et al. 



2004; 



Zhong et al. 



2003[ ). 



2001b 



Hahm and Liebei , 



On top of that, giant piezoresistance eff e ct ( 



2006 ) and enhanced thermoelectric performances (IBoukai et al. 



20081: 



He and Yang . 



Hochbaum et al 



20081 ) have recently been reported. The in t ereste d reader is encouraged to check some 



Patolsky and Lieber . 



2005 



reviews ( 


Kumar 




200'; 




Li et al 




2006; 


Lu and 


jiebei . 


Thelander et al. 




2006 




Wu et al. 


2008a: 


Xia et al. 


2003) 



2006 : 



In this paper we will review the theory of SiNWs. Clearly, we will make several references 
to experiments, whenever they support or challenge the theoretical predictions. Sometimes 
the comparisons are difficult to make, because SiNWs that are routinely grown range from 
50 to 200 nm, while those that can be efficiently studied within electronic structure methods 
are 2-3 nm thick, at most. Luckily, this gap is slowly narrowing and t hin SiNWs with 
diameters below 



2001a 



Cui et al 



nm have 



2001a 



Morales and Lieber 



1998 



jeen successfully 



2003 



De 



Wu et al. 



'adova et al. 



2004 



grown 



2008 



Zhong et al. 



by several groups ( 



Holmes et al. 



2000 



Coleman et al. 



Ma et al. 



2003; 



20051 ). The theoretical results that 



we discuss outline the most urgent problems that will have to be dealt with within the next 
generation of nanowires, those with characteristic sizes approaching the quantum limit. 

Although many of the features that we will discuss are common to other types of semicon- 
ducting nanowires, for the sake of clarity we will restrict to SiNWs throughout the paper. It 
should be at least pointed out, however, that in recent years tremendous progresses are be- 
ing made with compound semiconduc tors nanowires -mainly III-V nanowires- especially for 



what concerns photonics application (iBjork et al. 



20031). 



2002 



Dick et al. 



2004 



Thelander et al. 



A final remark concerns the computational methodologies. Although the main goal of 
the paper is giving a complete overview of the most important results that have been ob- 
tained within atomistic simulations, we will not enter in technical details, unless where 
it is necessary. Most of the results have been obtained within density functional theory 



(DFT), whose theoretical grounds are clearl y out of the scope of t 



reader can look at both t he original papers flHohenberg 



19651 ) , excellen t reviews (|Jones and Gunnarssonl . 



hensi ve books (iMartinl. 



ism ( 1 Colombo 



200 



atomic potentials (jJusto et al 



1989 



and Kohn 



lis work. The interested 
Kohn and ShanJ . 



1964 



Payne et al 



19921 ) and compre- 



2004r). Less frequently, we will refer 



Goringe et al 



1998 



1997 



Slater and Kostei 



Stillinger and Webei 



:o th e tight-binding formal- 



1954) or t o em pirical inter- 



1985 



Tersofa . ll989l ). 
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II. STRUCTURAL PROPERTIES 



A. Growth orientations and monocrystallinity 



The extra o rdina r y impact that the discovery of car bon nanotubes ( llijima . 



Oberhn et al. 



1976 



1991 



Radushkevich and Lukyanovichl . 



1952) had on condensed matter and 



nanoscience at first biased the research on Si quasi one- dimensional syst ems to the pur - 
suit of tubular structures. Hollow stru ctures resembling carbon n anotubes (iLi et all |2002| ). 



structures based on hoi 



ow elements (IMenon and Richter 



19991 ) or on fuUerene-like sys- 



tem (IMarsen and Sattlerl . Il999l) have been pr oposed. Although these -or other structures 



inspired by cluster assemble ( ISen et ali |2002| )- are stable within a total energy framework, 
they have not been observed experimentally to date. 



In the meanwhile Si nanotubes have been indeed successfully s ynthesized (ISha et al. 



2002 ). while their use for nanoelectronics still remains troublesome (jPerepichka and Rosei 



20061 ) . and things with Si nanowires turned out to be simpler than speculated. Convincing 
experimental evidence soon indic ated that SiNWs are rod-like structures constructed around 



a bu 



2003 



k Si single - crystalline core 



Wu et al. 



2004 



Holmes et al. 



Zhang et al 



200o[ ). 



2000 



Morales and Lieberl . 



1998 



Teo et all 



An important consequence of their single-cr ystal natur e is th at SiNWs grow along very 



Wu et al. 



(I2OO4J ) carried out an interesting 



well defined crystalline directions (see Fig. [2]). 
and extensive study of the growth orientations, showing a connection between the diameter 
and the favored crystal axis in the Au-catalyzed synthesis of SiNWs: the smallest-diameter 
nanowires grow primarily along the (110) direction, whereas larger nanowires favor the 
(111) direction; intermediate diameters, 10 to 20 nm, on the other hand, are dominated 
by (112) wires. Thermodyn amic models have been proposed to acc ount for this diameter- 



dependent growth direction (ISchmidt et al. 



2005 



Wang et al. 



2006al ) with consistent results 



in good agreement with the experiments, fixing the cross-over from (110) to (111) growth 
at 20-25 nm ((112) orientation was not considered in those studies). The stacking sequence 
preference that leads to (1 10) over (111) SiNWs at small diameters is also supported by 



first-principles calculations (lAkiyama et al 



20061 ). More recently, a continuum model that 



allows stu dying how growth begins an d evolves toward steady-state wire growth has been 
presented (jschwarz and Tersofi . 2009). The advantage of this approach is that complex 
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situations such as catalyst coarsening and interrupted growth can be easily handled. 

Ideally, nonetheless, one would like to be able to control the wire orientation at growth 
time. An important achievement in this sense was the demonstration th at the growth 



orientation can also 



2000 



3 e cont rolled externally by adjusting the growth pressure (IHolmes et al. 



Lugstein et al. 



20081 ). Alternatively, the use of different techniques can bias somehow 



the growth along certain crystal axis. For instance, the less common ox i de-as sisted growth 



method, generally yielding a broader diamete r distribut i on (|Wang et al. 



different orientations for ultra-thin SiNWs (ITeo et al. 



19981 ). might favor 



20031 ). Significantly, the thinnest 



SiNW reported to date (IMa et a/.l . l2003l ) was synthesized with this technique and was a 
(112) wire (see Fig.E}. 



B. Surface reconstructions in pristine nanowires 



The next major issue one has to face when studying the structure of a SiNW is the shape 
of its cross-section which, as we shall see briefly, is intimately related with the growth orienta- 
tion. Although one can pictorially imagine nanowires as cylindrical structures, clearly, when 
going down to the atomic-scale detail, this is not the structural arrangement that they as- 



sume -or even can assume. 



1984 



Zhdanov and Kasemi 



le an alogous problem in solids and small particles (jWang et al. 



19981 ) is elegantly solved by means of the Wulff construction or 



Wulff rule (iMarkd . Il994j ). which relates the equilibrium shape with the surface free energy 
of the facets involved. Solving the energy minimization problem min^s7s, where s is the 
nu mber of surface unit cells a nd 7^ the corresponding energy, leads to the optimum shape. 



Zhao and Yakobsoru (120031 ) have reexamined the use of Wulff construction within the 



determination of the equilibrium cross-sections of SiNWs. They showed that the conven- 
tional formulation of the Wulff criterion lacks of two important aspects: (i) in solids and 
smooth spherical particles the energy of the edges between facets is neglected compared to 
the surface contribution; (ii) the bulk is assumed already at its minimum and thus invariant. 
Hence, they propose the following generalization for the Wulff energy: 



(1) 



where they include the energy of matching adjacent facets E^, i.e. the energy of the edges, 
and the energy of the bulk E^, releasing the constraint on the innermost part of the wire 
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which can now change. 

They investigated different faceting arrangements for SiNWs grown along the (110) axis 
comparing them on the basis of Eq. [T] They found that the ground-state structure for SiNWs 
up to 5 nm is a pentagonal cross-section constructed joining five prisms cut out of a [110] Si 
plane [see Fig. Eld)]. This structure ha s seldom been detected experimentally [a remarkable 



observation by 



Takeguchi et al\ (l200ll ) is shown in Fig. [5](e)]), probably because it is not 



constructed around a bulk-core, which seems to be the favored situat ion at growth time 



Howev er, if one restricts to wires with a strictly bulk core the model of 



Zhao and Yakobson 



20031 ) correctly predic ts hexagonal over square cross-se ctions for (110) SiNWs, in agreement 



with the experiments (IMa et al. 



2003 



Wu et al 



The most important result of the work of 



20041) fsee Fig. 171). 



Zhao and YakobsonI (120031 ) is emphasizing the 



role of the edges and how the interplay between edges and surfaces play a key role in de- 
termining the reconstruction of Si one- dimensional structures. Befor e its formalization, this 
effect had been already pointed out by llsmail-Beigi and Ariad (119981 ) a few years earlier. In 
their work they considered a pristine SiNW grown along the (100) axis. This orientation fa- 
vors a square cross-section with jlOOj facets, an energetica lly cheaper solution than a square 
cross-section with {111} facets ( iRurali and Lorentd . l2005bl ). The abrupt match between the 
{100} facets results in an energetically expensive edge, a large value of Ee in Eq. [1], which 
can be reduced by forming smaller, transition {110} facets that allow a smoother match 
between the dominant {100} facets and partially release the stress accumulated at the edge 
(examples can be seen in Figs. 111(c), [6](a) andO). 

Unfortunately, a word of care should be spent concerning the above discussion. Down at 
the ultimate nanoscale limit it is delicate to give general rules and for extremely thin SiNWs 
counterexamples can be found to the general trends discussed pre viously. For instance, 
( I2OO6I ) showed that the faceting arrangement prop osed bv llsmail-Beigi and Arias 



Cao et al. 



2005 



(|l998l) for (100) wires and later followed by oth er authors (ILee and Rudd 



2007b 



Rurali and Lorente 



2005a 



Vo et al. 



Ruralil . 



20061 ) is favored only beyond a 1.7 nm diameter. 



whereas tiny SiNWs prefer s harp edges, i.e. the removal of the ed ges does not pay back. 



In the spirit of the work of 



Zhao et al. 



Justo et al. 



(|2007l ) carried out an interesting 



(12004J ). 

and systematic study of SiNWs grown along the (100), (110), ar id (112) crystal ax is, carrying 
out extensive calculations based on an interatomic potential (jJusto et all Il998l ). In order 
to elucidate the role of the different facets for the stability, for each growth orientation they 
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examined cross-sections bounded by different facet compositions. For the (100) SiNWs, 
for instance, they considered both all-{100} facets, all-{110} facets and three intermediate 
combinations. Proceeding in this way they were able to formulate a universal scaling law 
in terms of the wire perimeter, according to which the nanowire energy per atom always 
lies within two limiting energy lines, which are directly related to the character of the 
prevailing facet s. Interestirigly, in the limit of thick wires, the edge energy become negligible 



as suggested by 
perimeter. 



Zhao et al. 



(I2OO4J ) and the energy scales linearly with the inverse of the wire 



Silicon has a very rich phase diagram (IKaczmarski et all l2005l ) and many solid 
phases other than the diamond structure are known. Among them is the so-called 
clathrate phase that becomes stable at negative pressures. The stability of such a phase 
for quasi one- dimensional nanost r uctur e has been investigated by Ponomareva and co- 



workers (IPonomareva et al. 



20061 . I2OO5I ). They studied cage-like SiNWs carved out of a 



Si clathrate structure and compared them with both tetrahedral diamond-like and poly- 
crystalline SiNWs. Their results indicate that also in these nanostructures the tetrahedral 
structure is favored. Nevertheless, the difference in energy is rather small and it is suggested 
that clathrate based SiNWs might have better conductive properties. 



C. Passivated nanowires 



The study of the structure of pristine SiNWs has been a fertile ground where to start the 
theoretical research of these fascinating systems. However, quite soon it became clear that 
the wires grown experimentally have always passivated facets. Silicon form highly directional 
covalent bonds according to the know sp'^ tetrahedral pattern. Silicon atoms at the surface 
have dangling bonds (DB), unsaturated bonds that make the atom highly reactive and that 
induce strong reconstruction of the surfaces. Generally speaking, surface passivation consists 
of the termination of DBs on the surface with elements that assure their chemical stability. 
Hence the surface is chemically passive. 

Surface passivation in SiNWs mainly originates from two causes: (i) the growth of a 
thin layer of Si02 by thermal oxidation of silicon; (ii) presence of hydrogen in the growth 
environment during the synthesis or HF attack of the oxidized wires after growth, a process 
yielding removal of the Si02 layer and H passivation. Hydrogen passivation is rather simple 
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to model. If a sufficient amount of hydrogen is supplied the H atoms readily terminate 
each Si DB by forming a stable Si-H system. Passivation by oxidation is more complex. 
Thermal oxide is amorphous and then difficult to model at the nanoscale, because of the 
large amount of atoms required to describe the disordered phase. In the study of SiNWs, 
for most practical effects, hydrogen termination is a reasonable approximation to oxide 
passivation and this is the strategy adopted in most of the theoretical studies reviewed 
here. This approach is also justified by the fact that it is easy to remove the oxide layer 



after the growth a nd to induce H passiva t 



is often followed (IGuichard et al. 



2005 



2006 



Wang et aL i2008i : iWu et al. 



ion by simply etching it with H 



He and Yang . 



2006 



Ma et al. 



. Th i s procedure 



2003 



Ross et al. 



20041 ) in order to work with cleaner structures where 



the passivation rely on an individu al termination of the DBs, rather than a less controllable 
and more defective oxide coverage (IBaumer et al.\ . \2004 ). Furthermore, it has also proven to 
leave the morphology of the nanowire essentially intact, except for the removal of the oxide 
layer (jzhang et al . 2000), allowing inspection of the underlying atomic scale structure (see 
Fig. [3]). Yet, more attention is likely to be devoted in the near fu ture to the spec i fic na ture 



of Si02 passivation, beyond the simple models considered so far (lAvramov et all 120071 ). 



The passivation has a crucial effect on the electronic structure of the wires and it is essen- 
tial to provide the wires with predictable band gap widths and an invariably semiconducting 
character. We will discuss these topics in detail in Section IIIII 

Surface passivation has also an important effect on the structural arrangements of SiNWs. 
Besides preventing complex reconstructions, it also influences the structure of the sub-surface 
and innermost part of the wires. H-passivated SiNWs grown along different orientations 
have been found to maintain remarkably the bulk symmetry (see Fig. Hj), with negligible 
deviations of the Si-Si bond lengths; the deviation incre ases close to the surface, depending 
on the level of surface rearrangement (jVo et all |2006| ). The limit case in this sense are 
pristine wires, where the absence of passivation results in maio r surface rearrangemen ts and 
large deviations of the Si-Si bond length also in the wire core ( iKagimura et a/.l . l2005l ). 



An interesting pat 
has been proposed by 



1 to the determin ation of the structure of H-passivated (110) SiNWs 
(120061 ). Their optimization procedure is based on a genetic 



Chan et al. 



algorithm. With this method, in principle suitable for any other growth orientation, they 
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identified a pool of magic structures ^ for (110) wires. Although some of them have not 



been observed experimentally, thei r hexagoii a 



the STM image of the wire facet of | Ma et al. 

A more systematic app roach was 
to the one performed by lJusto et al. 



structure provided a good agreement with 



(2003L 



lowed by 



Zhang et al. 



( I2OO5I ) in a study analogous 



(120071 ) for pristine nanowires. They carried out a 



comprehensive study of the possible low-index facets in H-passivated SiNWs grown along 
the (100), (110), (111), and (112) axis. While many choices are possible for (100), (110), 
and (111) wires, they showed that there is only one low index configuration -with two {111} 
and two {110} facets- for (112) wires. It is suggested that this would ease the controlled 
growth with a predetermined cross-section and could have important consequences on the 
engineering of devices based on SiNWs. In Section IIII.B[ however, we will see that it has 
been recently suggested that the exact cross-section shape is less important than other 
parameters -such as the effective diameter and the surface-to-volume ratio- when it comes 
to determining the electronic properties of SiNWs. 

Another important aspect to consider in H-passivated S iNWs is t he su rface structure 



of the hydrogenated facets. This issue has been tackled by IVo et al\ ( 120061 ). where a sys- 
tematic study of the effects of varying the diameter and the growth direction has on the 
structure of the hydrogenated surfaces of (100) (110) and (111) SiNWs. In particular, they 
studied the relative stability of symmetric Sill2 dihydrides, canted SiH2 dihydrides and a 
(2x1) surface reconstruction (where first reconstruction is allowed and then passivation oc- 
cur), see Fig. El They deliberately chose simple, round cross-sections, as their scope was 
focusing on the atomic scale structure of the facet. Their wires were constructed selecting 
all the atoms falling inside a virtual cylinder placed in bulk silicon, in such a way that 
the facets approximated a circular cross-section. This procedure agrees with the smooth- 
ness prescription des cribed above, w hich -more importantly- seems also to be confirmed 
by the e xperiments (|Ma et alx . l2003l ). They found that, in agreement with bulk Si (100) 
surfaces ( iNorthrupl . Il99ll ). the canted dihydride surface is more stable than the symmetric 
dihydride structure, because canting allows a larger H-H separation. Additionally, faceting 
confers an increased stability to the canted dihydride surface, because at the facets' edges the 



Magic structures is used in this context to refer to distinct types of wi re configurat i ons w ith low formation 



energies that emerge as the number of atoms per length is increased (jChan et al\ . 120081 ) 



11 



SiH2 groups are free to rotate. R elief of the sur f ace st rain through bending as an additional 



mechanism has been explored by 



Zdetsis et al. 



mm . 



A perhaps more flagrant effect of the surface induced strain i s the fact that the axial 



Ng et al. 



(120071 ) reported 



lattice parameter of thin SiNWs is in general different from bulk Si. 
contraction along the wire axis for (100) (111) and (112) SiNWs, and elongation for (110) 
growth orientation. 



Other types o 



Nolan et al. 



surfac e passivation -inc ludi ng OH (lAradi et al. 



2007 



Ng et al 



2007; 



20071), NH2 flNolan et adl2007h . F flNg et adl2007h . or Br, CI, and I flLeu et al 



20061 )- have been considered. While changing the passivation has a limited effect the struc- 
tural properties of the nanowire, it can affect in a more significant way the electronic band 
structure. We will come back on this topic in Section IIII.B.3I 

As a conclusive remark one should notice that, despite the intensive research carried out to 
find the equilibrium shap es for the different growth orientations -pro posing structu r es tha t 
range from fuUerene-like (IMarsen and Sattled . Il999h to star-shaped ( ISorokin et all l2008h - 
in most cases the experimentally observed cross-sections of passivated SiNWs are decep- 
tively simple (see Fig. [3]and[71), whereas unpassivated SiNWs have never been reported. 
Furthermore, as we shall see in Section lllll although the cross-section shape has captured 
great attention and has been the object of many studies, in realistic, passivated wires the 
growth orientation and the average diameter turned out to have a more significant impact 
on the electronic properties of SiNWs. 



D. Mechanical properties of nanowires 



If one carves out of bulk Si a rod-shaped system like a nanowire, there is no apparent 
reason to expect an enhanced stiffness, while the larger surface-to-volume ratio is rather 
suspected to be detrimental. A simple way of understanding these effects is that there is a 
layer of material at the surface and edges whose mechanical properties differ from those of 
the bulk including different elastic moduli and eigenstrains. 



These intuitive ideas have been rigorously tested by iLee and Ruddl (l2007bl ) , by means 
of an exhaustive study of (100) SiNWs with increasing diameters. They calculated the 
Young's modulus, finding that it softens from the bulk value as the surface-to-volume ra- 
tio increases, going through a steep decrease around 2-2.5 nm diameter (see Fig. [8]). They 
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showed that the origin of this behavior is the compressive surface stress. To get a better 
insight into these atomic scale mechanisms the Young's modulus can be decomposed into a 
core (Si core atoms) and a surface contribution (Si surface atoms, H-H and Si-H systems). 
This decomposition allows highlighting the insensitivity to the facet ratio, as the contri- 
butions to the Young's modulus that are strongly facet dependent are very small. These 
first-principles re sults are in good agree ment with empirical atomistic potentials and contin- 



2007al ). unless for the smallest wires where these simplified 



uum techniques (ILee and Ruddl . 
approaches fail (see Fig. [8]). 

The Young's modulus, as many oth e r pro perties of ultra-thin SiNWs reviewed in this 
paper, is strongly anisotropic. iMa et al\ (120081 ) extended the study of the stiffness vs. diam- 
eter to wires grown along the (1 10), (111), and (112) o rientations. While their results are 
in good agreement with those of Lee and Rudd ( 2007b ) for (100) wires, they showed that 
wires of similar diameter, but with different orientations, differ considerably. In particular 
they found the highest values for (110) SiNWs, while (100) SiNWs giv e the lowest values. 

(I2OO8I ) where the Poisson 



Leu et al. 



These results are in good agreement with the work of 
ratio is also considered. 

We note that to calculate the Young's modulus a definition of the cross-sectional area 
must be assumed and it is not univocal. We will run into a similar problem concerning the 
definition of the wire diameter in Section IIII.B.2I when dealing with quantum confinement. 



Lee and Ruddl (|2007br) to ok the area bounded by the outermost atoms, i.e. the passivating 



H atoms, 



Ma et al\ ( 120081 ) made a similar choi ce, but excluding the H atoms. Aware of this 



degree of arbitrariness in the possible choices. 



Leu et al. 



20081) studied the variation of the 



calculated mechanical properties as a function of the uncertainty dr^ in the estimation of the 
radius tq. The error in the Young's modulus is 25ro/ro and goes to zero in the limit of large 
wires (rg 00). The Poisson ratio, on the other hand, is much more sensitive, because the 
error is — (5ro/ro + {5r — 6ro)/{r — tq), r being the radius at a strain e; while the first term 
vanishes for large radius, the other is always present and can be significant as (r — tq) is 
typically small. 

Bending has received co mparatively less attention, although a fe A V experimental measure 



ments have 



j een re ported (IHoffmann et al. 



Zheng et al. 



2006; 



Hsin et al. 



2008 



Tabib-Azar et al. 



2005 : 



20091 ). This is probably due to the difficulty of studying a bended nanowire 



within atomistic simulations that normally relies on periodic boundary conditions. The fab- 
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rication of ingenious mechanical structures with enhanced elastic p roperties suggests that 



20071 ). 



this could be a promising research direction (iSan Paulo et all 

Beyond elastic deformation, materials undergo non reversible, plastic deformation 
which directly precede f racture. This regime has been st udied for (100) pr istine 
SiNWs f Justo et al . 2007 ) and for (111) and cage-li k e SiNWs (Menon et al. I2O O4), us- 



ing two different interatomic potentials (jjusto et al. 



1998 



Stillinger and Weberl . Il985l ). 



At small strains e the stress increases linearly, as expected in the elastic regime, while 
at larger deformation the plastic behavior appears until the fracture occurs at e ~ 0.10, 
with a good agreement between the two different models. Experimentally, however, the 
fracture i s much delay e d wit h respect to theoretical predictions and the wire breaks at 



e ~ 0.25 (IKizuka et al. 



2OO5I ). It should be mentioned that both these theoretical studies 



considered pristine nanowires, while the wires in the experiments are coated with a thin 
layer of oxide, thus it is difficult to make a rigorous comparison. At the same time it is 
not surprising that a different surface treatment can produce a noticeable difference in the 
mechanic response, because it is just at the surface that the nanoscale signature emerges. 
Correlating st ructural defo r mation s and changes in the electronic properties is an active 



field of research (IRurali et al 



been investigated (IHong et al. 



2008af ) and the use of strain to enhance ca rrier mobility has 



2008 



Huang et al. 



2008 



Leu et al. 



20081 ). Furthermore, a 



giant piezoresistance effect -the application of a strain to a crystal that res ults in a change in 



the electrical resistance- has been reported recently ( IHe and Yang) . l2006l ). The underlying 



atomic scale mechanism is still poorly understood, however, and the attempts made so far 



proved to be somehow elusive (ICao et al. 



20071 : 



Rowel . l2008f ). 



III. ELECTRONIC PROPERTIES 



A. Pristine Nanowires 



The reason for the attention devoted to geometrical features such as the growth orienta- 
tion, the faceting arrangement, and the surface structure, is that they are crucial when it 
comes to the electronic properties of the nanowire. Clearly, the thinner is the nanowire, the 
more it is sensitive to the structure details, as in the limit of very large diameter -no matter 
which is its crystal axis or cross-section shape- its properties converge to those of bulk Si. 
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As we discussed in the previous section, pristine nanowires turned out to have a hmited 
relevance, at least to date, because experimentally grown SiNWs are always passivated. 
However, the study of bare, unpassivated wires is still interesting for two reasons: (a) it 
leads to the important conclusion that passivation is essential to obtain nanowires with 
predictable and easy to control electrical properties; (b) it sheds a light on some atomic 
scale mechanisms of high fundamental interest. 

An interesting example in this sense is the electronic structure of (100) SiNWs with {100} 
facets. While other facets, like the I 111) fac e ts, have an e l ectron ic structure similar to the 



corresponding infinite surface f Pande 



1981 



Rurali et al. 



20061 ), {100} facets can be very 

inown 



19791) 



different. In the Si(lOO) surface each surface atom has two DBs. The surface is 

J 

to reduce its energy by forming dimers, thus halving the number of DBs ( IChadil . 
The reconstruction of {100} facets follo ws the same pattern , but w ires dominated by such 
facets have been reported to be metallic. iRurali and Lorentd (j2005al ) showed that thin (100) 
SiNWs sustain two different reconstructions of the {100} face t that turn the wire metallic o r 



semimetallic, in agreement with what previously suggested by 



Ismail-Beigi and Ariad (1l998l ). 



The metallic behavior can be ascribed to a modified coordination of the jlOOj facet atoms 



leading to a distortion of the surface dimers, with respect to the Si (100) surface (IRurali et al. 



20061 ). The metallicity rapidly vanishes as the diameter is increased and the facets recover 



the coordination and the semiconducting electronic structure of the S i (100) surface 



The (100) wires with sharp corners studied by 



Cao et al 



(I2OO6I ) can be metallic too. 



Interestingly, the edge metallic states decay slower with the diameter co mpared to the facet 



Rurali and Lorente 



metalli c states. Consequently, wires thicker than those considered by 
(l2005al ). where edges were absent, can still be metallic. In both cases the metallic states 
are related directly or indirectly with the edges -in one case purely edge states, in the other 
dimer rows with an altered coordination near the edges; as the wire size increases the relative 
number of atoms at the edges decreases rapidly and the facet recovers the semiconducting 
character of the infinite surface. 

Besides the fundamental interest of these findings -in one case a metallicity driven by 
the finite size of the facet, in the other sustained by the edges- it is clear that such wires 
are not desirable for electronics application. On the one hand one wants to work with 
semiconducting systems; on the other hand, although some application can be envisaged for 
metallic SiNWs, e.g. interconnects, the metallicity should be much more robust, so that it 
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is not destroyed by small variations of the diameter and does not depend critically on the 
atomic scale structure of the wire. 

A comprehensive study of the surfac e reconstruct i on an d electronic structure of pris- 



Singh et al. 



(l2005l ). The cross-section chosen for 



tine (110) wires has been carried out by 
these wires i s such that they have {100} and {110} facets, at variance with the SiNWs of 



Rurali et al. 



(120061 ) which have {111} and {110}. This variation results in significant struc- 
tural differences, because of the compa ratively large r readjustment of {100} facets, which 
involve the formation of surface dimers (jChadil . Il979l ) and therefore a noticeable reconstruc- 
tion; in {110} and {111} facets, on the other hand, no new bonds are formed and the overall 
reorganization of surface atoms is moderate. These (110) SiNWs turned out to be indirect 
band gap semiconductors, with the states of the top of the valence band and the states of 
the bottom of the conduction band originating at different facets. Yet, it should be noted 



Fernandez- Serr a et al. 



that a m etallic reconstruction for (110) SiNWs have been reported by 
( l2006bl ). Once again, small variations of the atomic scale structure or of the cross-section 
can result in major changes in the electronic structure. 

Contrarily to what we shall discuss in Section IIII.B.2I concerning quantum confinement, 
here the thinner is the wire, the smaller is the effective band gap. This looks like a general 



featu r e of pristine, unpass ivated nanowires, where band gaps are smaller than in bulk ( 



2005 



Rurali et al. 



Ruralil . 



20061 ). In thin wires the surface-to- volume ratio is larger and surface 
states, which often lie in the gap, dominate the electronic structure and result in an effective 
narrowing of the energy gap. 



B. Passivated nanowires 



1. Band structure and hand gap anisotropy 



We have already mentioned a few times that surface passivation is required to obtain 
ultra-thin nanowires that are semiconducting and have a predictable and controllable band 
gap. Notwithstanding, the electronic structure of the nanowires still depends on the growth 



orientation, on the cross-section s^ 



anisotropic 


Leu e 


t al. 


2006 


, 2008 


Singh et al. 


2006; 


Vo et al. 


2006; 



lape and on t 



Ng et al. 



le diameter. The 



2007 



Niauet et a, 



Yan et al.\ . i2007i : iZhao et al 



j and gap is strongly 



2006 



Rurali et al. 



2007a 



20041 ) and, for wires of com- 
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parable diameters, it follows the ordering 



^(100>>^(lll>_^(112>>^(110)^ 



(2) 



with the orientation effect still sizeable up to 3 nm diameter (INg et ali l2007l ). The band 
gap of (112) w ires is of t h e same order of (1 1 1) wir es, though it has bee n reported to b e 



slightly larger (ILeu et al. 



2006 



Rurali et al. 



2007al ) or slightly smaller (INg et al. 



20071 ). 



This anisotropy has been qualitatively tracked back to the different geometric al structure 



of the wires in the (100), (111), and (110) directions ( iBruno et al. 



2005 



2007b|). While the 



(100) and (111) wires appear as a collection of small clusters connected along the axis, the 
(110) wires resemble a linear chain (see Fig. [TT] where the case of Ge nanowires is shown). 
Therefore, one expects that quantum confinement effects are larger in the (100) and (111) 
wires, because of their quasi zero-dimensionality, with respect to the (110) wires. 

Bulk Si has an indirect band gap, with the valence band maximum at the F point and the 
conduction minimum at about 85% along the F to X direction, and a phonon is required to 
conserve the momentum in any electronic transition. Remarkably, however, SiNWs grown 
along most of the crystallographic orientations have a direct band gap, meaning that the 
maximum of the valence band and the minimum of the conduction band occur at the same 



point in /s-space. This property has a.. 



materials for photonics applications (jCanham 



lowed to envisage the use of SiNWs a. s optically active 



1990 



Guichard et al. 



20061). 



In (100) SiNWs, the confinement plane contain four of the six equivalent conduction 
band valleys. These minima at ±?/ and ±2; are then projected onto F due to band folding, 
thus resulting in a direct band gap. When the axis is along a lower symmetry direction 
the confinement plane cannot contain four conduction band valleys and it will contain at 
most two. This is the case of (110) SiNWs. The minima at ±z are projected onto F again. 
Now, both the large and the small masses appear in the confinement plane, with the larger 
longitudinal mass being the relevant effective mass for describing the confinement effect in 
the cross-section plane. On the other hand, the four remaining minima will be projected to a 
point between F and the zone boundary Z, with the effective mass on the confinement plane 
being a value between the longitudinal and transverse masses. Therefore, the conduction- 
band edge at F is expecte d to have a smal ler upward shift induced by confinement and the 



band gap becomes direct (lYan et al. 



20o3). 



Although the projection along the (111) axis would lead to an indirect band gap, 
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the thinnest (11 1) SiNWs have a direct band gap (IRurah et ali l2007al : IVo et all 12006 : 



Zhao et al 



20041 ) ■ One should bear in mind that, besides the band folding arguments given 



above, the effective masses play an important role. In the quantum confinement regime 
(see below. Section IIII.B.2P the conduction band states are shifted upward, the smaller the 
diameter, the larger the shift. Nonetheless, the magnitude of this energy shift is different 
for each Appoint of the band structure and depends on the effective mass. In bulk Si the 
effective mass at F is heavier than at X or L. Hence, upon confinement, one expects the 
conduction band energy at X and L to increase more than at F. This simple considerations 
based on effective mass theory (EMT) describe well the transition from direct to indirect 
band gap experienced by (111) SiNWs that occur around 2 nm: as the diameter increases 
the quantum confinement effect vanishes (see Section [111.6.2^ . F and X and L points are 
not shifted and the gap remains indirect. Among the studied orientations (112 ) SiNWs are 



the only ones that 



2007 



20051 ). 



Huang et al. 



lave an indirect band gap also f or the thinnest diameters (lAradi et al 



2008 



Lu et al. 



2008a 



Ng et al 



2007 



Rurali et al. 



2007a; 



Scheel et al 



Although the band gap is highly anisotropic and, as we shall see in the next section, 
strongly dependent on the wire diameter, it is very inte r esting to observe that it is rather 
insensitive to the shape of the cross-section. iNg et al\ (120071 ) have studied the effect of 
the variation of the cross-section in thin (110) SiNW, generating 13 different cross-sections 
obtained by modifications of a reference 1 nm wire. They found that the band gap is 
practically constant and changes are within 0.09 eV. Later, it was demonstrated that wires 
of even utterly different cross-sectio ns can have the same band gap, provided that their 



surface-to-volume ratio is the same (lYao et al. 



20081 ). The effect of the surface-to- volume 



ratio on the band gap can be described by the universal expression 



E. 



gap 



(3) 



where E''^^^ is the gap of bulk Si, a is an adjustable parameter and S is the surface-to- volume 



m nm 
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2. Quantum confinement 



One of the most intriguing physical effect that arise in confined systems hke SiNWs is the 
so-called quantum confinement. Such a regime is conveniently descri bed through the particle- 



m-a- 



Dox model system in most quantum mechanics text books (IBransden and Joachainl . 



2000). The simplified situation considered is an infinite potential well where the motion 
of the particles is restricted to be in the direction of the confinement. As the motion of 
the particles is restricted, their kinetic energy increases and it is readily shown that the 
eigenstate energies are given by the following relation: 



(4) 



2m* (P 

where m* is the effective mass and d the width of the potential well. According to Eq. HI 
not only the energy levels, but also the spacing between them increases as the confinement 
becomes more pronounced, i.e. the smaller is d. Quantum confinement has a critical impact 
on semiconductors because it affects directly their most important electronic property: the 
energy band gap. 

Semiconducting nanowires provide a very good approximation of the model situation 
described above. Clearly, the potential well is not infinitely deep and realistic wire cross- 
sections like those described in Sections III.BI and III. CI are difficult to describe analytically, 
thus there is a need for a detailed electronic structure modeling. 

The first experimental proofs of quantu m co nfinement in nanostru c tured Si were reported 



in the pioneering works of 



Canhanj (jl990l ) and iLehmann and Goseld ( 119911 ) , where a simple 



electrochemical etching process was used to create crystalline Si nano structures with visi 



ble lu minescence at room temperature. As TEM images revealed later (jCuUis and Canhanj . 



199ll ). the etched structures consisted of rather disordered bundles of nanowires, t hough it is 



interesting to note that ordered structures like those speculated in the first place (jCanhaml . 



1990) have been recently proposed for the fabrication of ordered ar rays of quantum 



wires (Rurali et al 



Read et al 



2007br) and to achieve enhanced thermoelectric effect ( Lee et al. 



(jl992|) and 



Buda et al. 



20081 ). 



(I1992I ) performed DFT calculations of the band gap 
up shifts in perfect H-terminated SiNWs as a function of wire thickness, modeling porous 
Si (ICanhaml . ll990l ) with rectangular columns oriented along the (100) axis. Both these works 
showed that the fundamental gap is direct at the F point. This makes by itself the probability 
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of radiative recombination higher than in bulk Si, since no phonon is required in the electron- 
hole recombination process. Unfortunately, as it is well-known, standard local and semi-local 
implementations of DFT fail to account quantitatively for the band gap of semiconductors 
and one must resort to self-energy corrections to the Kohn-Sham gap to obtain a good 



agreeme nt with the exper i ment al val ues. Yet, t 



correct (IWilliamson et al. 



20021 ) and 



Read et al. 



l e tren ds are expected to be qualitatively 
( I1992I ) reported band gap upshift of up to 



2 eV for wires of ~ 12 A diameter. They also shown that a generalization of Eq. Hlgives a good 
description of the quantum confinement for wires wider than 23 A, whereas t hinner wires 



Buda et al. 



show significant deviations from this idealized EMT picture. In such a range 
(119921 ) showed that with the more realistic DFT potential the band gap scales as the inverse 
of the diameter (i, rather than as predicted by particle-in-a-box arguments where in- 
finitely hard walls ar e assumed. Subsequent studies that the intere s ted reader could look at 
i nclud e the works by 



([1993|L 



and 



O 



mo et al. 



Xia and Chang 



Ossicini et al\ (119971 ) 



I993I ) 



1992), 



Sanders and Ghana (11992 



Yeh et al. 



mm 



Saitta et al. 



19961 ) 



ivbertsen and Needels 



Xia and Cheahl (|l997f ). 



A first step toward a qua ntitative evaluation of SiNW band gaps in the quantum con 



finement regime was given by 



Delley and Steigmeierl (Il995l ). including a constant self-energy 



correction independent on the size. Namely, they increased all their calculated band gaps 
by 0.6 eV, the self-energy correction for bulk Si. They also shown that EMT can predict 
with great accuracy the band gap of relatively thin nanowires, provided that the potential 
well is not assumed to be infinite. 

The self-energy correction to the Local Density Approximation (LDA) or the Gener- 
alized Gradient Approximation (GGA) band -gap, howev e r, is e xpected to depend o n the 



wire diam e ter an d on the growth orientation. 



Yan et al. 



Zhao et al. 



mm 



Bruno et al. 



(l2007af ). and 



(120071 ) carried o ut calculations within the many-body perturbation method based 
on the GW approximation (lAryasetiawan and Gunnarssonl . ll998l ) for (100), (110), and (111) 
SiNWs. They showed that the self-energy is indeed anisotropic and is larger for thinner wires. 



The dependence of the band gap on the wire diameter can be described as 



E,a, = E^l' + C X (l/ci)^ 



(5) 



where E^J^^ is the calculated band gap of bulk silicon, d is the effective diameter of the wires, 
while C and a are fitting parameters. This formula is derived within a simple particle-in- 
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a-box effective mass approximation, where a = 2 when barrier height is infinite. The GW 
results can be fitted to this formula, yielding values of a ranging from 0.9 to 1.1, much lower 
than those expected within EMT and depending on the growth ori entation (see Fig. [TUD 



so that the band gap and the dielec tric response are anisotropic (IBruneval et al 



Bruno et al. 



2007a 



Zhao et al. 



2005 : 



200J). 



Although GW is in principle the best suited methodology to calculate the band gap in 
semiconductor systems, it s uffers from the seri o us inconvenience of a cons iderable co r nputa - 



tional load. In the works of 



Zhao et al. 



mm . 



Bruno et al. 



(12007ah and 



Yan et al. 



(120071 ). 



for instance, only relatively small SiNWs can be calculated directly and the band gaps of 
larger, more realistic wires are obtained by numerically fitting the available data to Eq. O 
Furthermore, an alternative to many-body GW calculations is mandatory when it comes to 
calculate doping levels, a task that requires large computational cells. In the remainder of 
this section we discuss two possible approaches. 

A successful way to improve the DFT band gaps consists in using hybrid-functionals for 
the exchange-correlation energy, where a certain amount of exact Hartree-Fock exchange 
is mixed to conventional LDA/GGA functionals. The amount of Hartree-Fock exchange 
(typically 12-15%) is chosen to reproduce some parameters of the bulk system (the band 
gap, among them), rather than being based on solid theoretical grounds. Hence, strictly 
speaking, one cannot claim to solve the electronic structure from first-principles. In such 
a theoretical fr amework th e band gap of SiNWs with diameter s up to 3 nm can be calcu- 



lated directly ( lAradi et al. 



2007 



Ng et al. 



2007 



Rurali et al. 



2007al ). These results are 



important becaus e they allow dir ect comparison with the only experimental measurements 



available to date (IMa et al. 



20031 ). A direct comparison of the experimental data with GW 



calculations is not possible for two reasons: the diameters of the wires grown experimentally 
are larger than those that could be simulated and most of the available measurements are for 
for (112) SiNWs, whose larger primitive cell precludes GW calculations even for the thinner 
wires. 

Alternatively, the band structure of nanowires can be calculated with a semiempirical 



tight-binding method, where the sel 



ment of the image charge effects (INiquet et al. 



-energy is ob t ained within a simpler semiclassical treat- 



2006[ ). 



This is a very powerful method 



because, due its reduced computational load, it allows calculating SiNWs with diameters up 
to 10 nm with good accuracy. As we will see in some more detail below (Section IIII.C. 31) . a 
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great advantage of this method is that it allows dealing with different dielectric surr oundings 



which is very important in systems with abrupt dielectric interfaces like nanowires (iLi et al. 



20081). 



So, how should the band gap of SiNWs be calculated? The accurate calculation of band 
gap is one of the most challenging problems in semiconductor theoretical physics, so it 
is not surprising that it is not easy to answer this question. GW calculations provide in 
principle the most accurate estimations. However, they are restricted to very thin SiNWs. 
Semiempirical tight-binding, on the other hand, is a very attractive choice for larger wires, 
especially for those diameters where quantum confinement become small and gap broadening 
is dominated by dielectric mismatch effects (jPereira et all 120091 ). Hybrid- functional DFT 
calculations are an interesting compromise for those wires that are too large for GW (too 
many atoms in the primitive cell) and too small for tight-binding (relaxation effects cannot 
be neglected and the use of a parametrization obtained for bulk Si could be questioned). 
It is difficult to asses on the accuracy of each of these methods, since the experimental 
measurement of the energy gap of SiNWs is extremely challenging, and only the data of 



Ma et al. 



(120031 ) are available to date. More experimental results are needed to clarify this 



important point. 



3. Surface chemistry 



We have seen above -and will see again below- that many properties of SiNWs are deter- 
mined by their large surface-to-volume rat io. Hence, i t is na t ural that most of the exc iting 
physics takes place at the wires' surface (IKobayashil . 



2004 



Zhong and Stocks! . l2006l ) . In 



Section IIII.AI we have seen, for instance, that wires bounded by facets derived by semicon- 
ducting surfaces can exhibit surface metallicity. Passivated nanowires are more predictable, 
in this sense, and it is just because they are always semiconducting that they are expected 
to play an important role in the next generation of electronic devices. Yet, the surface has 
a relevant role that merits some considerations. 



An important case is that of chemical sensors whe re the adsorption of a mo 
yields measurable variations of the electrical conductance ( iBlase and Fernandez-Serral . 



ecule 



Cui et al. 



2008 : 



2001bl ). Upon adsorption, the molecular orbitals can hybridize with the wire 



states, resulting in a sizeable modifications of its electronic structure. Whether the ef- 
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fectiveness of this process depe nds on the facet where adsorption takes place is the topic 



addressed by 



Leao et al. 



(120071 ) ■ They studied the sensitivity of different facets of a (110) 
SiNW, showing the existence of a specific relation between the way surface atoms are bonded 
to core atoms and the relative contribution of these surface atoms to band edge states. These 
observations are important concerning the optimal design of those chemical sensors where 
the adsorption of a molecule directly modify the wire transmission. A broader class of sen- 
sors, however, seems to work on simpler basis. The dipole induced by molecule adsorption 
can act as a gate voltage, opening or closing the conductive channel in a field-effect transistor 
set-up. 

In the quantum confinement regime the band gap width depend critically on the di- 
ameter. The possibility of controlling the band gap width is tremendously attractive for 
optoelectronics applications: not only SiNWs can have a direct band gap, which per se 
increases the optical efficiency, but its width can in principle be tuned. It is not difficult 
to imagine, however, that controlling the wire diameter with tolerances within 1-3 nm is a 
more than challenging task. A simpler route to band gap tuning is controlling the chemical 
composition and the coverage density of the wire surface. Halogens such as CI, Br, and I 
can be used as surface passivation agents instead of H and, while not altering th e semicon- 



ducti ng character of the wires, they result in a significant shrinking of band gap (ILeu et al. 



20061 ). The strongest reduction of the band gap is provided by I, followed by Br and CI, in 



the opposite order of the bonding strength of these species and SiNWs. Interestingly, the 
surface coverage is a further degree of freedom and one can span all the band gap values 
between a H- and halogen-passivated wire by varying the H:halogen ratio. Also, increasing 
the halogen surface concentration the band edge states, concentrated in the wire core in 
presence of H-passivation, progressively spread to the surface^ 



Analogous results have been reported for OH and NH2 (lAradi et al. 



2007 



Nolan et al. 



20071 ). It should be noted that the passivating species do not contribute significantly to 



the states close to the band edges, so that the reduction of the gap is not caused by the 
introduction of additional bands. It rather comes from the hybridization of the valence band 
states with the frontier orbitals of the different passivating functional groups that cause a 
significant band gap reduction relative to H-passivated wires. 

These results indicate that the band gap width in SiNWs can be tailored not only by 
controlling the wire diameter, but also by an appropriate choice of the surface termination. 
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C. Doped and defective nanowires 



Semiconductors are privileged materials for electronics applications because their resis- 
tivities can be can be varied by design with great control ^. Equally important, they can be 
designed to conduct one of two types of carriers: electrons and holes. These two features are 
the core of device design, which relies on the interaction of adjacent semiconductors with 
different densities of electrons and holes. The most efficient way to control the carrier den- 
sity is doping the semiconductor, that i s incorporating substitutional impurities (dopants) 



in the lattice (IMuUer and Kamind . Il986l ) . 



In the simplest model, a substitutional defect with minor relaxation effects forms four 
bonds with the neighbor atoms in the Si crystal. For a group-V element, such as P or As, 
the fifth valence electron is not covalently bonded to near neighbors and it is only weakly 
bonded by the excess positive charge of the impurity nucleus. Hence, a small amount of 
energy is required to break this weak interaction and this electron is free to wander about 
the crystal and contribute to conduction. These impurities are called donors because they 
donate an electron; analogous arguments apply to group-Ill elements which are acceptors. 

There are at least two reasons that make the physics of impurities in nanowires different 
with respect to bulk systems: (i) the lattice sites are no longer equivalent in the direction 
of confinement and (ii) in the quantum confined regime all the eigenstates, including those 
associated to defects, are shifted in energy with important consequences on the dopant 
activation. Below, we discuss in some detail these and other topics relevant to dopant 
efficiency. We conclude the section with a generalization of the formation energy for defects 
in nanowires. 



1. Surface segregation, surface traps and dopant aggregation 

In bulk Si all the lattice sites are equivalent. In a nanowire this is true only for the 
axial direction, because the lateral confinement breaks the translational symmetry. In other 
words, given one arbitrary Si atom in the nanowire, as one moves along the wire axis he finds 



^ Conventionally, materials with resistivities less than about 10 ^ cm are consid ered conductors and 
mate rials with resistivities greater than about 10^ cm are considered insulators i Muller and Kamins , 



19861) . 
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an infinite number of equivalent atoms, whereas as one moves along the radial direction he 
finds a finite number of nonequivalent atoms. Therefore, the substitutional defects at these 
non-equivalent sites will have in gene ral different formation e nergies and doping levels. 



This problem was first tackled by 



Fernandez- Serra et al. 



fl2006bh . who studied B and P 



substitutional in (100) and (110) SiNWs. They revealed a tendency to surface segregation of 
these impurities, which means that substitution close to the surface is energetically cheaper 
than substitution in the innermost part of the wire. The effect is especially pronounced 
in presence of DB defects, so that the dopant impurities are effectively trapped by these 
surface defects. Most importantly, the dopant-DB complex is electrically inactive, reducing 
the carrier concentration at room temperature. 

The segregation energy of P is significantly larger than that of B. This means that for 
the same impurity concentration, a much larger fraction of P atoms will be captured and 
neutralized by surface tr aps, result i ng in a lower conduct ance. This is in agreement with 



the experimental results (jCui et al. 



2000 



Yu et al. 



2000) where, for similar doping levels. 



B-doped SiNWs present a lower resistance than P -doped ones. Similar studies have bee n 



carried out for (110) SiNWs of different diameters fiLeao et al. 



2008 



Peelaers et al. 



20061). 



They showed that B and P prefer to sit at edge or near edge sites (the most external lattice 
sites with all-Si nearest neig hbors) , depending on the surface facet and on the atomic surface 



structure (ILeao et al. 



20081). 



In perfectly passivated wires the surface segregation originates from a simple relaxation 
effect. At these edge and near edge positions it is easier to release the strain induced by 
the substitution, whereas in the center of the wire the allowed relaxation is smaller due to 
the constraint of the surrounding Si lattice. In surface defective wires, on the other hand, 
the presence of DBs greatly enhances the tendency to surface segregation, with the impurity 
atoms moving at surface sites. Here the driving force to surface segregation, yielding a much 
more sizeable effect, is an electronic effect: the formation of a stable dopant-DB complex. 

The surface trapping of dopants has a dramatic impact for two reasons: (a) the impurities 
trapped at the surface are deep-state defects and are electrically inactive, thus not contribut- 
ing to the carrier concentration at room temperature; (b) due to their large surface-to- volume 
ratio and considering a typical dopant concentration of 5x10^^ cm~^ and an estimate of of 



10^^ cm ^ interface DB defects, one finds that for wires of 



always enough DBs to trap all the dopants (jPeelaers et al 



ess than 4 nm diameter there are 



20061 ). Luckily, the difference in 
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formation energies among surface and core substi t utiori al sites has been shown to decrease 
rapidly as the diameters grow larger (ILeao et all l2008l ). Hence, there is a twofold benefit 
in enlarging the wire diameter: on the one hand, the surface-to-volume ratio decreases, and 
so does the density of DBs with respect to the dopants concentration; on the other hand, 
the trapping efficiency of these reduced density of surface defects is lower, as the forma- 
tion energi e s of c ore substitutionals and dopant-DB complexes become of the same order. 



Leao et al. 



(120081 ) estimate that the populations of core and surface dopants will be similar 



for wires with diameters of 3 nm or more. 



The importance of surface impurities has been further highlighted by iDurgun et al. 
( l2007al ). They considered various impurities such as Al, Ga, C, Si, Ge, N, P, As, Te, 
and Pt, focusing on adatoms configurations, rather than on substitutionals. They found 
that the energetically most favorable adsorption site of the six considered depends on the 
group of the Periodic Table that the impurity belongs to. All the configurations studied, 
however, give rise to deep state in the gap and are not viable choices as active dopants. 

Another source of dopant deactivation is the formation of electrically inactive dopant 
complexes. Two nearest neighbor dopants can form a bound state, so that the weakly 
bonded electrons that contributed to the conductance are now participating in the dopant- 
dopant bond. Yet, this energy gain has to compete with the energy cost that results from 
the strain accumulated around the dopant-pair defect. This str ain is more easily r eleased in 
nanowires than in bulk, due their large surface-to- volume ratio 



Moon et al. 



(120081) reported 



a high stability of P pairs, which increases as the wire diameter is reduced. Interestingly, this 
is not the case of B. When two B atoms occupy nearest neighbor sites the lattice undergoes 
a significant relaxation, the B impurities move far apart and assume a planar threefold 
coordinated configuration. This is possible because, unlike P, B can present either sp^ or sp^ 
hybridization. Again, ^^type doping seems to be more robust than n-type doping, at least 
as far as one considers B and P. A similar mechanism leads to mutual passivation when both 
a B and a P impurity are present at the same time. B esides the obvious compensation of 



having an n- and a jo-type dopant, 



aggregation at the wire surface. Als o, the use o : 



Peelaers et al\ (j2006l ) showed that these impurities favor 
B/P co-doping has been proposed to reduce 



in a controllable way the band gap (jlori et al. 



20081 ). 



26 



2. Quantum confinement 



In Section IIII.B.2I we have seen how one of the most important quantities of a semi- 
conductor, its band gap, depends critically on the wire diameter in the regime of quantum 
confinement. This is a very important parameter, because it determines the amount of 
carriers -intrinsic carriers- that can be thermally excited from the valence band to the con- 
duction band. Intrinsic carriers, however, are not very important at typical device operation 
temperatures and the conduction is dominated by extrinsic carriers, those carriers that are 
thermally excited from a dopant level. Hence, dopants to be such must be very shallow, 
meaning that the impurity electronic states have to be only a few meV from the band edge. 

Now, if the band gap broadens as an effect of quantum confinement, what happens to 
dopant levels? In a purely effective mass picture they will be shifted, like any other state, 
becoming deeper, i.e. the upshift of a donor level will be less than of the conduction band 
edge. Clearly, this fact has dramatic consequences on the dopant efficiency. Namely, a 
dopant impurity which is known to be very shallow in bulk Si, becomes deeper as the 
diameter shrinks down, and it will not eventually be usable to dope ultra-thin SiNWs. At 
which diameter does this happen? 

From EMT, one can deduce the effective Bohr-radius of the ground state as ~ (e/m*)ao, 
where Oq is the Bohr-radius of the isolated hydrogen atom. This results in about 2.2 nm 
(thus a 4.4 nm diameter) for P. A crude estimate of the extension of the wave function is 
taking twice this diameter, thus ~ 9 nm. Yet, EMT neglects relaxation effect which can be 
important in very thin SiNWs and the dopant levels should be calculated directly. 

The trend of th e ionization e n ergies v s diameter can b e qualitatively obtained from 



DFT calculations flDurgun et all l2007a] : iLeao et all 120081 ). As we discussed in Sec 



tion nil. B. 21 however, the local and semi-local approximations commonly used for the 
exchange- correlation energy severely underestimate the band gap and likewise the gap states 
and the related ionization energies. While the best suited approach for correct band gap 
calculations was the GW methodology (see Section [III.B.2I) . it does not seem a viable solu- 
tion for defective systems. Due to the need of simulating isolated impurities, computational 
supercells have to be large enough to allow neglecting the interaction of a defect with its 
periodic images. This implies a large number of atoms which is normally beyond the current 
computational capabilities of GW based codes. 
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Hybrid functionals -where a certain amount of exact Hartree-Fock exchange is mixed to 
con ventional LDA/G GA functionals- provide accurate estimations of defect states in bulk 



SiNWs f lRurah et al 



Si (IDeak et all 12005 ) and have also been used to calculate P donors in (110) and (111) 



20091 ). As expected, these calculations yield ed ionization energies 



that are deeper than the values obtained by DFT (ILeao et alx . |2008| ). though the difference 
decreases for larger wires. Remarkably, P behaves as an EMT dopant down to diameters 
of 1.5 nm, its wave function highly localized, whereas it breaks down for wires of 1.0 nm 
diameter. For such small wires the wave function is qualitatively different: it significantly 
interacts with the surface and cannot be described as quasi-one dimensional confined EMT 
state. 



3. Dielectric confinement 



The estimation of the ionization e nergies of dopants in nano wires has also been tackled 



efficiently at the tight-binding level (jPiarra et al 



2008 



20071 ). We have already seen in 



Section III1.B.2I that thi s approach can be complementary to the calculations at GW and 
hybrid-functional level flNiquet et a/.l . l2006l ). Its quantitative reliability could be questioned 
for ultra-thin wires (diameters smaller than 2 nm), because this model neglects relaxation 
effects, important for such wires, and relies on a parametrization obtained for bulk Si. On 
the other hand, it is the best alternative to deal with larger wires where confinement still 
produces sizeable effects. A remarkable feature of this approach is the flexibility with which 
the screening properties of the surrounding dielectric medium can be manipulated, allowing 
to study in detail the so-called dielectric confinement. 

The Coulomb potential of an impurity gives rise to a bound state in the energy gap. 
In bulk Si this potential is strongly screened (e^ = 11-3), the Bohr radius is large, and 
the ionization energies amount to a few meV, so that the impurities are ionized at room 
temperature. The screening of a nucleus charge +e leaves a total charge +e/e at the impurity 
site, whereas the remaining charge +e(l — 1/e) is repelled at infinity. In a one dimensional 
system the screening properties are different. The charge +e(l — 1/e) must be repelled at 
the surface of the nanowire, leading to an extra term in the potential. 

The physics of screening in one dimensional systems is straightforwardly incorporated in 



28 



the tight-binding Hamiltonian of iDiarra et al\ (120081 . 120071 ) 



H = Hq + Uimp + ^ 



(6) 



where Ho is the Hamiltonian of the undoped wire, Uimp = ^V{r,rQ) is the potential at r 
of an impurity at tq, and ^ is the self-energy potential, which accounts for the interaction 
between the carrier and the surface polarizati on charges in d uced by its own presence. 



On the basis of this tight-binding model 



Diarra et al. 



(120071 ) showed that the ioniza- 



tion energies of typical donors are significantly deeper than in bulk, even for large wires 
{d > 10 nm) where the effects of quantum confinement are weak. This effect is due to the 
interaction between the electron and the surface polarization charge +e(l — 1/e). These 
results (i) indicate that dielectric confinement can be stronger than quantum confinement 
and that donor levels deepen more than how much the band gap broaden; (ii) the dielec- 
tric mismatch can be used to vary the ionization energies. In particular, a metallic or 
high-permittivity surrounding gate, present in realistic applications, is expected to reduce 
significantly the ionizat ion energies . Thes e predictions have been recently supported by the 
experimental results of iBjork et all ( 120091 ). 



4- Metallic impurities 



Although much of the attention has been devoted t o dopants so 



are becoming increasingly irnportant for the reliability (IBailly et al. 



ar, metal impurities 



2008 



Hannon et al. 



20061: 



den Hertog et al. 



2008 



Oh et al. 



20081 ) and functionalization of SiNWs. Transition 



metals have attracted some inte rest because of their possible use in de s igning na r iosca. 



dilute magnetic semiconductors (IDurgun et al. 



2008 



Giorgi et al. 



2008 



Xu et al 



2008). 



Room-temperature ferromagnetism in SiNWs has indeed been reported (iWu et ail 120071 ). 
although the annealing conditions to stabilize the magnetization are very critical. The 
ferromagnetic coupling of Mn impurities is confirmed by electronic structure calculations. 
At variance with typical d opants th e se ini purities do not segregate to the surface, at least 
in absence of surface DBs f iXu et all 120081 ). and they favor aggregation. This tendency has 
been reported to be important to stabilize the magnetism, because only when the Mn-Mn 
distance i s below a certain c utoff ferromagnetic coupling is favored over antiferromagnetic 



coupling (iGiorgi et al. 



20081). 
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It is interesting to observe that specific TM impurities can drive the nanowire to a half- 



metaUic grou nd state. A hal: 



is insulating (Ide Groot et al. 



- metal is a system where one spin is metallic, whereas the other 
19831 ) ■ Such systems have the greatest interest for spintronics 
applications, because naturally all the conduction electrons belong to the same spin and the 
spin polarization P is maximum ^. I n the framewo r k of a comprehensive analysis of surface 
adsorption of TM atoms in SiNWs, iDurgun et al\ (l2007bl ) discovered that wires decorated 
with Co and Cr can be ferromagnetic half-metals. Upon adsorption of a Co or Cr atom 
at the surface, the spin degeneracy is lifted and, while the bands of majority spin continue 
being semiconducting, two minority spin bands made of hybridized TM-3d and Si-3p states 
cross the Fermi level driving the wire to half-metallicity. There is a sizeable charge transfer 
from the TM atom to the wire and the high values of the binding energies indicate strong 
bonds, which is important to prevent uncontrolled clustering that would be detrimental to 
the magnetic ordering. The half-metallicity is obtained for huge coverages (one impurity for 
primitive cell) typical of dilute magnetic semiconductors. As the coverage is reduced, the 
gap of the minority spins starts closing and the system is no longer half metallic {P < 1). 
However, the total spin polarization remains very high and close to its maximum permitted 
value. 



5. Formation energy 



The energetic cost of creating a defect is given by the formation energy. The forma- 
tion energy is the main computational quantity describing the stability and energetics of 
a defect in a host material ari d it is essential to determine impurity e quilib rium concen- 



t r at i ons (IZhang and Northrup 



ties (IFahey et al. 



19911), solubilities (IVan de Walle et ali Il993h or diffusivi 



19891 ). In bulk sy stems the formation e n ergy can be calculated according 



to the well established formalism of 



Zhang and Northrupl (jl99ll ): 



(7) 



where E^^^ is the total energy of the defective system, Ui is the number of atoms belonging 
to species i with chemical potential /ij, q is the net charge of the system, is the energy of 



^ The spin polarization is normally defined as P = {N'^ ~ N^)/{N^ + iV^), where A^T^i ^j-g the densities of 
states at the Fermi level. 
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the top of the valence band of the clean host and /ig is the chemical potential for electrons; 
the sum runs over all the species present in the system. However, the extension of Eq. [7] to 
one dimensional systems presents some subtleties. In particular, the chemical potential of 
the host species is ill-defined. 

To see why, let us focus on the formation of a Si vacancy and let us play the movie of 
the wire growth: Si atoms start precipitating from the supersaturated Au-Si droplet and 
contribute to the nucleation of the wire; once in a while one Si atom does not fill the proper 
lattice site and a vacancy is formed. When calculating the formation energy one has to 
estimate the contribution to the total energy of this atom, the one that left the vacant 
site and was incorporated in the wire somewhere else. In bulk this is easy, because all the 
lattice sites are equivalent and then each atom contributes equally to the total energy. In 
a nanowire, on the other hand, this quantity is not well defined, because it depends where 
th e misplaced atom is thought to end up, as the lattice sites are non-equivalent. 

Rurali and Cartoixal ( 2009 ) have proposed a way to circumvent this problem that also 



deals with the passivating agents, if present. They showed that rewriting the equations for 
the formation of defects at a time, being N the number of Si atoms in the primitive cell, 
leads to the definition of an effective chemical potential of the wire primitive cell. In this 
way, whenever the formation of a defect involves the addition/removal of an atom of the 
host species, e.g. vacancies, self-interstitials, substitutionals, it is transferred to/from the 
correct reservoir -the wire itself- and only easily computable quantities are involved. 

A further problem arises when dealing with charged defects. In a periodic boundary con- 
dition formalism point charges result in an infinite electrostatic energy. This inconvenience 



can be obviated by using a compensating jellium background. Th e errors in the tota. 



are often corrected a posteriori by means of the Madelung energy (IMakov and Paynd . Il995l ) , 



energy 



though other methods have been proposed. Again, this correction has to be properly general- 
ized to be used in nanowires. In particular, in solids the Madelung correction is scaled by the 
value of the (isotropic) macroscopic dielectric constant of the host material. In nanowires, 
on the other hand, a dielectric tensor e will be needed for the correct description of the in - 
teraction between the different instances of the charged defect (iRurali and Cartoixal . 120091 ) . 



Notice that the value of the dielectric tensor will depend on the ratio between the axial lat 



tice p arameter, the lateral vacuum buffer and the chosen values of the axial (IHamel et al. 



20081 ) and transverse components of the e tensor and therefore cannot be looked up in tables. 
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IV. TRANSPORT PROPERTIES 



The study of electron and heat transport is one of the most rapidly growing research 
field in nanowires. The reason is twofold: on the one hand transport measurements often 
are the most direct and simp 



Chen et al. 



2008 



Cui et al. 



est way to test the theoretical pre d ictions (iBiork et al 



2000 



Li et al. 



2003 



Sellier et al. 



2006 



Yu et al. 



2009 



2000); on 



the other hand, the behavior can be much different from bulk Si and can be exploited for 
enhanced performances in applications, whereas other times it can be detrimental. 



A. Electron transport 

1. Surface roughness disorder 

An important cause of the degradation of the electrical conductance in SiNW-based de- 



vices is the scattering occurring at the su r 



roughness (ILuisier et al. 



2007 



Wang et al. 



ace in presence of surface defects or surface 



20051 ). This is not unexpected, since we have 



seen that many properties are ruled by the large surface-to-volume ratio of SiNWs. Con- 



pering flKodambaka et al 



2006; 


Wang et al. 


2006c: 


„^ — „ — , 
Wu et al. 


2008b) 



2008b!) or fancier saw-tooth 



faceting (IRoss et alx . l2005l ) are often reported. 

The effects o f surfa ce roughness on electron transport have first been addressed by 
(120071 ). They modeled the surface disorder adding with probability 1/2 



Svizhenko et al. 



one monolayer at each facet of a given unit cell. Since the position of these surface bumps 
is uncorrected, they obtained a white-noise roughness. The surface roughness originates 
strong irregularities in the density of states along the wire axis, which in turn causes re- 
flection of carriers and a strong reduction of the conductance. When these effects sum up 



m very 



ong wires 



;he d isorder quickly drives the transport into the Anderson localization 
19581 ). 

A compl ementary approach t o the description of t he su rface roughness disorder has been 



regime (lAndersonl . 



followed by 



Persson et al. 



mm and 



Lherbier et al 



(120081 ). They modeled the roughness as 



random fluctuations 5r of the wire radius around its average value ro, through a Lorentzian 
autocorrelation function, obtaining a correlated disorder. They showed that the backscatter- 
ing strongly depends on the nanowire orientation, the anisotropy coming from the differences 
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in the underlying band structure. In particular, electrons are less sensitive to surface rough- 



ness in (110) SiNWs, whereas holes are better transmitted in (111) SiNWs ( jPersson et al. 



20081 ). Also, as the disorder correlation length -roughly the length scale of the diameter 
fluctuations- increases, the lowest-lying states of the conduction band get trapped into the 
largest sections of the wire ^. The modified extent of the electron wave function affect many 
key quantities for transport, such as the mean free path and the localization length. Inter- 
estingly, the room temperature mobility of electrons and holes seems rather insensitive to 
short length scale fluctuations, as well to very long length scale fluctuations, a case in which 
the surface experienced by the carriers is locally smooth. 

2. Single-impurity scattering 

Besides surface disorder, the other main critical source of scattering is the presence of 
impurities. Surface scattering has a stronger impact on the transport in SiNWs than in 
bulk Si because of the much larger surface-to-volume ratio. The case of impurity scattering 
seems different, since it should solely depend on the impurity density and should affect in a 
similar way bulk Si and SiNWs. This is not the case, though. So, where is the catch? 

With the reduction of the wire size below 10 nm, the impurity cross-sections become of 
the same order of the wire characteristic dimension and can result in total backscattering. In 
the semiclassical picture used to study transport in bulk materials impurities are point-like 
centers that scatter randomly the incoming carriers. This chaotic process slows down the 
carrier flow and results in a reduction of the conductance. The quantum picture in a thin 
one-dimensional medium is slightly different: impurities have to go through to a scattering 
potential that often extends throughout most of the wire cross-section and, following with 
the semiclassical analogy, the trajectories of the carriers are not simply deviated, but they 
can be entirely backscattered. 

Impurity is a fairly generic denomination when referred to semiconductors. In fact, it 
refers to both undesired defects, by-products of an imperfect growth, and to dopants, which 
are purposely introduced to provide the material with tailor-made electric features. Clearly, 
this case is the most challenging: dopants increase the carrier density at device operation 



^ A similar mechanism has been reported recently for selectively strained nanowires ( Wu et al\ . 12009 ). 
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temperature, but at the same time might induce a s ignificant scattering which leads to a drop 



in the conductance. 



Fernandez- Serr a et al. 



(l2006al ) studied the resistance associated with a 



substitutional P impurity in the wire core, at the wire surface, and with a DB+impurity, a 



comp 



ex whose importance was discussed in a previous work of theirs (IFernandez-Serra et al. 



2006bl ). Resonant backscattering -a strong reduction of the conductance in correspondence 
to impurity related bound states- is the main signature of substitutional impurities, though 
P in the core or at the surface yield different results. On the other hand, DB+impurity 
complexes are transparent to the incoming carriers and the transport is ballistic. 

Therefore, donor impurities such as P either segregate to the surface where they are likely 
to form an electrically inactive complex with a DB or they stay in the wire core where they 
produce a strong backscattering, particularly at certain resonant energies. In both cases the 



current is reduced. 



3. Multiple-impurity scattering 

The calculations of Fernandez-Serra and co-workers opened the field of dopant scatter- 
ing in SiNWs, but they have two limitations: (i) they study the scattering properties of 
an individual impurity, while in realistic SiNWs the wire resistance results from multiple 
scattering events; (ii) impurities can be ionized, the typical situation for dopants, and the 
proper charge state must be taken into account in the conductance calculation. 

The first of these two issues has been tackled by comparing the conductance evaluated 
directly in long wires, with a certain distribution of imp urities, with the predic tions that 



can be made on the basis of single-dopant calculations (IMarkussen et ali 120071 ). This is 
a challenging task, because to make such a comparison on equal footing the long wire 
too has to be treated within a first-principle formalism, which involves an extraordinary 
computational load. This difficulty can be circumvented thanks to an ingenious method that 
allows constructing the Hamiltonian of the long wire a ssemb ling building blocks obtained 



from the single-dopant calculations (IMarkussen et al\ . [200^. In this way the electronic 



structure problem has not to be solved directly in the long wire. 

The surprising result is that the properties of long, realistic wires -such as mean free 
path, resistance vs length- can be entirely predicted from single-impurity conductances. So, 
the resistance of a wire with an arbitrary distribution of impurities is obtained by classically 
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adding the resistances of each individual scatterer according to Ohm's law: 



{R{L,E)) = R,{E) + {Rs{E))L/l 



(8) 



where {Rs{E)) is the average resistance that can be evaluated from the single-dopant cal- 
culations, Rc{E) is the contact resistance, L is the length of the wire and I the average 
separation betwee n dopants. Interestingly, a similar approa ch has been proposed also for 



phonon transport (IMarkussen et al. 



200S 



Savic et al. 



20081). 



This method allows easy comparisons of the conductance of wires with different dis- 
tributions of def e cts. T he case of P substitutionals, for instance, has been addressed by 



Markussen et al. 



(l2008bl ). where a uniform radial distribution was compared to a mainly 



surface distr ibution, in accordance with the previously reported indications of P surface 



segregation (IFernandez-Serra et al. 



2mm 



4- Charged impurity scattering 

Addressing charged impurities poses well-known problems related to the use of periodic 
boundary conditions. Large supercells, out of the current capability of first-principles meth- 
ods, are needed to allow the correct screening of the electrostatic potential of the impurity. 
The conductance depends more critically than other quantities on this incomplete screening. 

Such systems have been dealt with within an approximate method that combine first 



princip le methods with finite element calculations of the electrostatic potential (iRurali et al. 
2008bl ). The idea is very simple. If a charged dopant is approximated with a point charge, its 
electrostatic potential can be obtained in a very cheap way with a finite element calculation. 
Far from the impurity this is a reasonable approximation -a P"*" impurity gives rise to 
essentially the same Coulomb potential of a As^ impurity- and the agreement with a self- 
consistent electronic structure calculation is indeed very good. Close to the impurity, on 
the other hand, quantum electronic structure accounts properly for the different chemical 
nature of different impurities, a task not accomplished by a finite element model. 

The part of the potential that converges slowly with the cell size is the long-range Coulomb 
potential. So, here comes the simple idea: the local potential around the impurity is cal- 
culated at the first-principle level in a large, but tractable computational cell. Then the 
long-range tails are calculated within a finite element model (which it is known to yield the 
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same result than a first-principle calculation). In this way one can in principle engineer 
arbitrary boundary conditions, as the long range part of the potential can be obtained at a 
negligible computational cost. 

This approximated approach allows explicitly addressing the calculation of the conduc- 
tance associated to impurities with different charge states and comparing majority carriers 
(electrons in a n-type wire, holes in p-tjpe wire) vs minority carriers (electrons in a p-type 
wire, holes in n-type wire). The results are indeed utterly different and it i s shown that 



in suff iciently thin wires minority carriers transmission is entirely suppressed ( iRurali et al. 



2008bl ). What happens is that in the case of minority carriers the dopant constitutes an 



effective barrier in the potential landscape. When the energy of the electron is less than that 
of the barrier height, it must tunnel through the potential and the transmission is therefore 
exponentially suppressed. 



B. Heat transport 



Recently, there has been a lot of excitement around the thermal conductive properties of 
SiNWs. Surprisingly, this excitement stems from the poor thermal conductance of SiNWs. 
The reason is that the use of SiNWs as materials wit h enhanced thermoelectric properties 



has b een demonstrated independently by two groups ( iBoukai et al. 

2mm. 



2008 



Hochbaum et al. 



While in some devices one wants to get the heat away as efficiently as possible and a 
high thermal conduct ance is desirable, i n thermoelectrics one wants a thermal conductance 
as small as possible ( Cahill et al. . I2OO3 ). It has been suggested that at the nanoscale the 
thermoelectric efficiency could be increased with respect to bulk materials, since the electrical 



mobilities are expected to be higher, whil e the sur: 



decrease the lattice thermal conductance (jVo et al. 



' ace sc attering of the phonons should 



20081 ). Recent experimental results 



have reported thermal conductivities of ~ 1.6 W m~^ K"^, two orders of magnitude lower 
than the value for bulk Si (150 W m^^ K^^ at room temperature). 

Heat is transmitted by phonons, the vibrations of the crystal lattice. Calculating the 
phonon modes of a SiNW with a first-principles me thod is a demanding task and can be 
done only for the thinnest wires (jPeelaers et a/.l . 120091 ). Fortunately, this is not too serious an 
inconvenience and the phonon band structure can be calculated with a great level of accuracy 
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within simple empirical interatomic potentials. It has been sho wn that phon on dispersions 



calculated with DFT and with the bo nd-order Tersoff potentia l (ITersofi 



conductances in excellent agreement (IMarkussen et al. 



19891 ) yield thermal 



2008al ) . It is indeed on the Stillinger- 



Weber potentia l that the first atomis tic calculations of the thermal conductance of SiNWs 



were based on (IVolz and Chen 



19991). 



The decrease of the thermal conductance with the reduction of the diameter comes from 
the interplay of two 



tra (Adu et ali 



"actors: (a) phonon confinement, i.e. the change in the phonon spec- 
20051 ) and (b) the increase of the inelastic phonon scattering at the sur- 
face. The dep endence of the p h onon dispersion on the wire size has been studied for 
(100) wires by iWang and WangI (120071 ), who showed that the thermal conductance de- 



creases as the wire diameter is reduced. 



Ponomareva et al. 



( 120071 ) obtained similar results 



for (111) wires, although the t herni al conductance steeply increases again for diameters be- 



low 2 nm (IPonomareva et al. 



J, 120071 ). a direct signature of phonon confinement. Namely, as 
the diameter is reduced the lowest frequency excited mode is severely affected by the con- 
finement; this long wavelength mode dominates the low frequency spectrum and, carrying a 
larger amount of energy, determines the enhanced thermal conductance at small diameters ^. 

It should be bear in mind, however, that the effects of phonon confinement are normally 
studied in nanowires with an ideal structure, whereas in the recent reports of the enhanced 
therm oelectric figures of S i NWs surface corrugation seems to play an important role. Re- 
cently, iDonadio and Gallil (120091) showed that the computed thermal conductance strongly 
depends on the surface structure, whereas it can be insensitive to variations of the diameter 
in the size range investigated (rf < 4 nm). Phonon confinement do not necessarily lead to 
low values of the thermal conductance, which in some cases can even increase as a function 
of size, d ue to the presence o f long wavelength phonons with very long mean free paths, like 



shown by 



Ponomareva et al. 



(120071 ). 

It is useful to co nsider these studies toge ther, because of their complementary method- 
ological approaches. IWang and WangI (120071 ) calculate the phonon dispersion with the usual 



procedure consisting in displacing each atom along ±x, iy, and ±2; to obtain the dynam- 



A similar feature is also present in the (100) wires of 1 Wang and Wang (|2007t ). where for temperatures 
lower than 50 K the thermal conductance of a 1.54 nm wire is larger than the thermal conductance of a 
4.61 n m wire. However, i t is dif ficult to be more quantitative at this respect, due to the scale of the plots 
used in 



Wang and Wantd (|2007f ) 
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ical matrix by finite differences. Then they calculate the thermal conductance with non 
equilibrium Green's functions. In this way they rely on the harmonic approximation and 
thereby neglect any phonon-phonon scattering, an approach vali d in the low to mid temper- 
ature limit. Similar calculations for thicker wires, d > 35 nm (iMingo 



2003 



Mingo et al. 



20031 ). where phonon conf inement effects are unimportant, y i elded an e xcellent agreement 



with e xperimental results (ILi et aUl2003l ). 



Ponomareva et al. 



mm and 



Donadio and GaUi 



( 120091 ) ■ on the other hand, calculate the thermal conductance within a molecular dynamics 
simulation from Fourier's law = — adT/ dz, where Jy is t he heat flow along the wire axis 
z and dT/dz is the thermal gradient ( jSchelling et all 120021 ). Within this scheme, one does 
not calculate explicitly the full phonon dispersion and anharmonic effects are automatically 
included in the simulation. The anharmonic forces are increasingly important at higher 
temperatures, since the atomic displacements get bigger. This means that phonon-phonon 
scattering becomes more and more important at higher temperatures and dominates over 
the effects of including more conducting channels. The drawback is that classical molecular 
dynamics is reliable only above the Debye temperature of the material (645 K for Si) where 
quantum effects in the ionic dynamics can be neglected and below it the results must be 
interpreted with some care. 

It should not come as a surprise that the thermal conductance is also anisotropic, like 
many other important quantities that we have discussed throughout the paper. At low 
energy the phonon disper sion features four acousti c bran ches, one dilatational, one torsional, 
and two fiexural modes (IThonhauser and Mahanl . 1200J). The torsional mode is related to 
rotational invariance around the wire axis, and it is similar for all the orientations. At higher 
energies, however, one can notice that the bands in (110) SiNWs have a larger slope, i.e. 
larger velocities, than the other orientations which feature mostly flat bands (see Fig. [T2l) . 
Hence, at a given energy there are more bands in the (110) wire than in the (100) or (111) and 
consequently one expects a larger thermal conductance. Heat conductance is indeed strongly 
anisotropic. Up to ~ 20 K, where the phonon dispersion is dominated by the acoustic modes. 



the thermal conductance is independent on the growth orientat ion, but then (1 



stand out, with an up to twofold increase in their conductance (IMarkussen et al. 



0) SiNW s 



2008a'). 
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V. CONCLUSIONS 



In this paper we have reviewed the major advances in the theoretical study of the struc- 
tural, electronic and transport properties of silicon nanowires. While the geometry and the 
electronic structure of nanowires are relatively well-understood, many open questions remain 
on the possibility of effectively dope ultra-thin nanowires and on many of the atomic scale 
mechanisms ruling electrical current and heat transport. 

Silicon nanowires are rod-like system constructed around a single-crystalline core. The 
most important consequence of their monocrystallinity is that they grow along well-defined 
crystallographic orientations, and at sufficiently small diameters a strong anisotropy of most 
of their properties emerge: the band gap, the Young's modulus, the electrical conductance 
or the specific heat, just to name some, are different for wires grown along different ori- 
entations. The cross-sections are intimately related with the growth orientations -given a 
growth orientation only certain sets of bounding facets are allowed-, although their impact 
on the electronic properties of the wires seems limited, while other magnitudes such as the 
diameter or the surface-to- volume ratio have a greater influence. The band gap can be di- 
rect, opening the way to the use of Si in photonics, and can be tuned by varying the wire 
diameter, choosing the growth orientation or controlling the surface passivation. Extrinsic 
carrier conduction seems to be highly problematic and whether is feasible or not for ultra- 
thin nanowires is not clear yet. The reason is that dopant efficiency is bedeviled by multiple 
factors: surface segregation and clustering with consequent neutralization, deepening of the 
doping level due to dielectric and quantum confinement. 

Many promising applications have already been demonstrated. Although the nanowires 
used in these applications are smaller than any device that can be fabricated with lithogra- 
phy based techniques, they are still larger than those studied theoretically, where quantum 
effects leave their clear signature. Whereas it is clear that nanowires will play an impor- 
tant role in the next generation of electronic devices, it is difficult to say if the use of such 
ultra-thin wires will be practical. Many of the properties of these extremely thin nanowires 
pose severe technological challenges, but at the same time represent extraordinary oppor- 
tunities. Anisotropic band gaps that critically depend on the wire diameter are apparently 
incompatible with any standardized technological process, to give an example. However, 
once the growth orientation and wire thickness can be controlled with great precision this 
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would open up the possibility of band gap engineering, which would be extremely attractive 
for optoelectronics applications. 

Joints efforts in theory and experiments hold the key to nanowires' future. 
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FIG. 1 (Color online) (a) Cartoon and (b) experimental realization of a ZnO nanowire-based field- 
effect transistor with an all-around (or surrounding) gate. The channel length is 200 nm. From 



Ng et al 
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FIG. 2 (Color online) Transmission electron microscopy (TEM) image of a single crystalline SiNW 
grown along (a) the (111) and (b) the (112) axis. The high resolution TEM micrograph of the 
crystalline core shows clearly the Si(lll) and Si(224) planes, respectively, together with the Fourier 
transform of the image. In case of the (111) SiNW a magnifie d view of the s i dewal ls of the wire 



show the presence of Au-Si particles about 7 nm in size. From 



20081) 
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FIG. 3 (Color online) (a) Constant-current scanning tunneling microscope (STM) image of the 
{111} facet of a SiNW grown along the (112) direction; (b) Schematic view of S1H3 on Si(lll) 
viewed along the (111) direction; (c) Constant-current STM image of the {100} facet of a SiNW 
grown along the (110) direction; (d) Schematic view of the dihydride phase on Si (001). Red and 
large blue circles represent the H atoms and Si atoms, respectively. Small blue circles correspond 
to Si atoms on the underlying layers. The (112) wi re in pane l (a) , w ith a diameter of 1.3 nm, is 



the thinnest SiNW reported to date. Adapted from 
from A A AS. 



Ma et al 



((20031). Reprinted with permission 
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FIG. 5 Hexagonal (110) wire with four {111} and two 1 100) (a) unreconstructed and (b) recon- 
structed facets, (c) Energy of different types of wires (see lZhao and YakobsonI (|2003l )) as a function 
of their diameter d. The most stable structure for d < 6 nm is [solid pentagons in panel (c)] is 



sh own in panel (d) 



by 



Takeguchi et al. 



'ane 



(d) shows a hig h resolution TEM linage of a pei itagonal nanowire 



(|200lh . Adapted from 



Zhao and YakobsonI (j2003l ) and 



Takeguchi et al. 



grow n 



(|200lh . 
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FIG. 6 (Color online) (a) Cross-section view of 3 nm SiNWs grown along three different directions 
(100), (110), and (111)- (b) Side view of three different sur face struc t ures: in the last configuration 
the surface first reconstructs and then is passivated. From 



Vo et al 



(12002) 
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FIG. 7 (Color online) (a) TEM images of 3.8 nm SiNWs grown along the (110) direction, (c) high 



resolution TEM cross-sectional image, and equilibrium shapes for the (b) N W and the 



cross-sections predicted by Wulff construction. The scale bars are 5 nm. From 



Wu et al. 
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FIG. 8 (Color online) Young's modulus calculated within DFT as a function of wire size. For 
comparison values of continuum formula are also plotted. The solid curve E = Ej^f^ — C/w, 
where w is the width of the wire and C = 66.11 GPa/nm, is the best fit to a pure surface area 
to volume size dependence. (Insets) Cross sections of some of the SiNWs studied, where each Si 



atom is colored corresponding to its transverse rel axation in A. T 
(b) 2.05, (c) 2.80, and (d) 3.92 nm. Adapted from 



Lee and Ruddl (I2007bl 1. 



le wid ths of wires are (a) 1.49, 
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k k k k 

FIG. 9 (Color online) Band structures of (100) (110) (111) and (112) SiNWs with a diameter of 
~ 3.0 nm (cross-sections in the insets). The arrows indicate the fundamental band gap which is 
direct for (100) and (110) SiNWs and indirect for (111) and (112) SiNWs. As discussed in the text 
the band gap of (111 ) SiNWs becomes direct when the diameter is reduced below 2 nm. Adapted 



from 



Ng et al. 



2003). 
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FIG. 10 (Color online) Quasi-particle GW gaps for (100) (circles), (111) (squares), and (110) 
(diamonds) SiNWs as a function of wire size com pared with experimental results (triangles) from 



scanning tunneling spectroscopy (iMa et al 



20031 ) ■ T he gray regio n 



gaps from (110) (bottom) to (100) (top) wires. From iBruno et al. 



repres ents the LDA electronic 



(|2007al l. 
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FIG. 11 (Color online) Geometrical structures of 0.4 nm Ge nanowires along the (110) (top), (111) 
(middle), and (100) (bottom) directions shown from the side (left) and from the top (right). Large 



spheres repres ent Ge atoms ; smal 
Adapted from 



Bruno et al. 



(120051). 



spheres are hydrogen atoms used to saturate the dangling bonds. 
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FIG. 12 Phonon band structures calculated within the Tersoff bond-order potential of 2 nm diam- 
eter SiNWs grown along the (a) (100) (b) (110) and (c) (111) axis. The phonon wave vectors, q, 
are ah in the respective wire directions and are sho wn in units of the recipr ocal unit ceh lengths, 



with aioo = 5.4, ono = 3.8, and am = 9.4 A. From iMarkussen et al 



(|2008al ). 
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